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Abstract. The Dyson-Schwinger equations of Landau-gauge Yang-Mills theory for the gluon and ghost 
propagators are investigated. Numerical results are obtained within a truncation scheme which has proven 
to be successful at vanishing temperature. For temperatures up to 250 MeV we find only minor quantitative 
changes in the infrared behaviour of the gluon and ghost propagators. The effective action calculated from 
these propagators is temperature-independent within the numerical uncertainty. 

PACS. ll.10.Wx Finite-temperature field theory - 12.38.Aw General Properties of QCD - 14.70.Dj Gluons 



1 Introduction 

It is by now well established that for increasing tempera- 
ture QCD undergoes a phase transition from a confining to 
a deconfined phase. Despite the phenomenological success 
of QCD, the understanding of the phenomenon of con- 
finement is still far from being satisfactory. The same is 
true for the deconfined phase: the picture of a quark-gluon 
plasma implies certain properties of quarks and gluons at 
high temperatures that have not yet been convincingly 
demonstrated. 

This state of affairs is due to the fact that both, con- 
finement and deconfinement, are not describable by per- 
turbation theory. Thus, investigations of this and related 
topics have to rely on non-perturbative methods. Lattice 
Monte-Carlo simulations of the Euclidean path integral 
provide the most direct approach and have provided con- 
vincing evidence that the pure SU(3) Yang-Mills theory 
undergoes a phase transition at a critical temperature 
T c « 270 MeV, see e.g. [Q. The confining properties as 
measured by the (temporal) Wilson-loop change across 
this transition which indicates that quarks are, at least 
partially, deconfined in the high-temperature phase. The 
fate of gluon confinement, on the other hand, is much less 
clear. 

At vanishing temperature, gluon confinement has been 
related to the infrared behaviour of the gluon and ghost 
propagators in covariant gauges, for reviews see e.g. (2IEI 
or for compact presentations of this issue 0). These prop- 
agators have been studied in lattice calculations, see 
and references therein, employing Renormalisation Group 
techniques [5] and Dyson-Schwinger equations (DSEs), 
see [7] and references therein. The results of these different 



methods agree qualitatively and quantitatively: The gluon 
propagator is infrared suppressed and the ghost propa- 
gator is infrared enhanced. From the latter property the 
Kugo-Ojima confinement criterion [5j follows. Further- 
more, it is equivalent to Zwanziger's boundary condition 
on the Gribov horizon [TJJ (see also |TT|) which guarantees 
that only field configurations within the first Gribov hori- 
zon | 12) contribute. In this approach, the occurrence of 
gluon confinement is thus related to the Gribov problem: 
due to the long-range nature of Faddeev-Popov ghosts not 
only Gribov copies are avoided but also the long-range 
propagation of transverse gluons is inhibited. It is there- 
fore interesting to learn how such a picture changes with 
temperature. 

Here, we will present the solution of coupled DSEs for 
the gluon and ghost propagators at non-vanishing tem- 
perature T^O. We focus on approaching the phase tran- 
sition from the low-temperature confining phase. An in- 
vestigation studying the high-temperature limit has been 
presented elsewhere . A solution on a space-time torus 
will be given because the use of a compact manifold pro- 
vides a natural infrared regulator (comparable to the one 
present in lattice calculations). 

This paper is organised as follows: first, we briefly out- 
line the finite-temperature framework for the DSEs. For 
the sake of making the presentation self-contained we dis- 
cuss in section[3]the truncation scheme employed. We then 
present the numerical results for the propagators. These 
are then used to calculate the effective action |14U15) in 
our truncation scheme. In the last section we conclude. 
Technical details are deferred to two appendices. 
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2 Gluon and Ghost Dyson-Schwinger 
Equations at Non-vanishing Temperatures 



Throughout this investigation the imaginary-time formal- 
ism, see e.g. |16U17j . is employed. This implies using peri- 
odic boundary conditions in "time" not only for the gluon 
fields but also for the ghost fields (despite their Grass- 
mannian nature). In four-momentum space this leads to 
a sum over discrete Matsubara frequencies oj n — 2nnT 
rather than a continuous integral: 



d A q 



(27T) 



T 



E 



d 3 q 

(2^)3 



(1) 



involving the three-momentum q. Here, we have chosen to 
study the system in the rest frame of the heat bath thus 
working in a formalism where Poincare covariance is not 
manifest. 

At non-vanishing temperature the full gluon propa- 
gator acquires a more complicated tensor structure. In 
Landau-gauge, which we use throughout this manuscript, 
the gluon propagator is purely transverse with respect to 
the gluon four-momentum which, at zero temperature, im- 
plies that it can be described by one scalar function. At 
non- vanishing temperatures two independent tensor struc- 
tures exist, one longitudinal and one transverse to the heat 
bath. The decomposition for the rest frame of the heat 
bath |16U17| is explicitely given by 



5 a 



(P Tt ,u{k)Z m (k ,\k\) 
+P L „ v (k)Z (k ,\k\)) 



Prij(k) = Si 



kj^ k j 



Pt on — Pt m — Pt m — 



P L ^{k) = P^(k) - P Tlw {k), = <W - 

hj = 1,2,3; n, v = 1, 2, 3, 4, 



kfj, k v 



(2) 



(3) 



with k 2 = fcg + k 2 . The ghost propagator is a Lorentz 
scalar and does not acquire further structures: 



D%(k) 



gab 

-^G(k ,\k\). 



(4) 



Another major change in the renormalised DSEs as com- 
pared to the ones for zero temperature is the substitution 
CJ. The equations for the ghost and gluon propagator, 
Da(k) and D^ v {k) , respectively, read: 



D-?;(k) = Z 3 (5^ 



— 1 

+ 00 



g 2 N c Z{T ]T I 



d 3 q 



n— — oo 



(2tt) 3 

iq ll D G {p)D G (q)G v {p, q) 

n— — 00 v ' 

r $ a (k, -p, q)D afi {q)D p<J (p)r l3au {-q,p, -k) 
+ -.., (6) 

with p = k + q. Contributions involving the four-point 
vertex are not explicitly given here. The colour structure 
is also suppressed, rfg denotes the tree-level 3-gluon- 
vertex, r a p 1 the full three-gluon-vertex function and G v 
the full ghost-gluon vertex. Here, the following renormal- 
isation constants appear: Z3 for the gluon wave function, 
Z 3 for the ghost wave function, Z\ for the 3-gluon vertex 
r a pj and Z\ for the ghost-gluon vertex G v . In Landau- 
gauge one has Z\ = 1 |18j . This non-renormalisation of 
the ghost-gluon vertex is crucial for the success of the em- 
ployed truncation scheme [T5ll2T)| . Following [2H a MOM 
renormalisation scheme will be used. Thus the renormal- 
isation constants depend on the renormalisation scale fj, 
and the ultraviolet cutoff, for details see below. 

At vanishing temperature, the DSEs on a torus provide 
a good tool to study finite volume and periodic boundary 
conditions effects and allow to compare to results of con- 
tinuum methods and the lattice. To extend this approach 
to finite temperature one can directly apply the proce- 
dure described in j22] by allowing for a different size of 
the torus in time-direction. The corresponding length is 
determined by the temperature, L — [3 — 1/T (feg = 1) 
with Lq <C L, L being the size of the torus in the spatial 
directions. Thus we substitute 



T 



E 

n— — ck 



d 3 q 
(2n) 



£ E E 



L 



(7) 



n=-oo j a ,j 2 ,33 



In order to solve the equations numerically, one needs an 
ultraviolet momentum cutoff A. As we use the same renor- 
malization procedure employed at T = this should be an 
0(4) invariant cutoff also at finite temperature. The inte- 
grals are replaced by the Matsubara sums, and the system 
of equations is solved self-consistently for a finite number 
of sampling points of the dressing functions. 

In the limit of infinite spatial volume, the effective tem- 
perature on the torus is 



D G \k) 



-Z n k 2 



+ ^cZ lT + f 

n— — 00 

ikfiD^k - q)G u (k, q)D G (q) 



(5) 



T = l/L . 



(8) 



For small tori this equality can be affected by finite size 
effects which will be discussed in the appendix B, confirm- 
ing JHJ to be a good estimate. 
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3 The truncation scheme 

Since the DSEs form an infinite system of hierarchically 
coupled integral equations, one has to truncate the sys- 
tem to make it tractable. We employ a truncation scheme 
which has been tested in various ways at zero tempera- 
ture |22| . The consequences of such truncations have been 
studied extensively in the literature, see ref. for a re- 
view. In the ultraviolet the truncation is fixed by requiring 
consistency with perturbation theory. In the infrared very 
likely only the behavior of the ghost-gluon vertex is rel- 
evant jH)ll20ll23| . A bare ghost-gluon vertex as employed 
here is in agreement with recent studies using DSEs [23 
and lattice methods [2E|. Thus the truncation is proba- 
ble even exact in the infrared, as has been argued using 
stochastic quantization |llj . Thus only at mid- momenta 
truncation dependent effects are relevant, but are strongly 
constrained. The results in the vacuum are furthermore in 
very good agreement with lattice calculations |26|. 

In this scheme only the propagators of the ghost and 
the gluon fields are determined self-consistently, while full 
vertex functions either have to be constructed or are re- 
placed by bare ones. As already mentioned, the gauge pa- 
rameter is fixed to the Landau-gauge value. The Landau- 
gauge is a fixed point of the renormalisation group and 
thus the gauge parameter is not renormalised. Further- 
more, contributions from the four-gluon vertex are ne- 
glected; this is justified for two reasons. In the first place, 
the tadpole diagram only contributes a scale-free, diver- 
gent constant (in Landau-gauge) in the perturbative regime, 
which drops out by renormalisation. Secondly, it can be ar- 
gued that the two-loop diagrams are sub-leading in the in- 
frared, if the 3- and 4-gluon vertices do not acquire highly 
singular dressing. Thus we are left with the 3-gluon vertex 
and the ghost-gluon vertex which are not determined by 
the propagator equations. Since the ghost-gluon vertex is 
not ultraviolet divergent in Landau-gauge, replacing the 
full vertex by the tree-level one preserves the one-loop 
anomalous dimensions of the dressing functions. In the 
truncation scheme used, we simply set 



and 5 = —9/44 being the anomalous dimension of the 
ghost propagator. In the vertex ansatz (|1(J|) we have cho- 
sen the linear combination of eq. I|12|) for the two-gluon 
dressing function, because it corresponds closest to the 
zero-temperature dressing of the gluon propagator. The 
vertex construction (|10() . taken over from vanishing tem- 
perature studies, should also be sufficient at finite tem- 
peratures, since it mainly affects the ultraviolet regime 
which is nearly temperature independent up to very high 
temperatures. 

The truncation described causes spurious longitudinal 
contributions to the gluon polarization. These spurious 
terms are quadratically divergent. It is possible to iso- 
late the part without quadratic divergences by contracting 
with the traceless Brown-Pennington projector |27| . 



(9) 



For the 3-gluon vertex the following construction pro- 
vides the correct one-loop behaviour in the ultraviolet 
for the propagators without perturbing their infrared be- 
haviour |H] 

r pua (k,p,q) = F(k,p,q)rW(k,p,q), (10) 
1 G(k)^ 2 ^ 



F(k,p,q) = 



Z x {ij?,A 2 ) Z(fc)(!+ 35 ) 
G(p)(- 2 - 6 ^ 



with 



Z(k) := l -Z Q {k) + 2 -Z m (k) , 

r pll(k,P,q) = -i{k-p)p5 liv - i{p-q)^ vp 
- i(q - tyvS^p , 



(11) 



(12) 



(13) 



K(k) 



i kpk u 



(14) 



which projects out terms proportional to 5^. However, 
this projector interferes with the infrared analysis. In or- 
der to take care of this problem at zero temperature, the 
gluon equation is contracted with the transverse projector 



V(k) = v - ^ 



(15) 



and a quadratically divergent tadpole- like tensor structure 
is subtracted from the gluon loop in the gluon DSE. It is 
given by [Hj 



Q 



5 1 



4 k 2 q- 



(16) 



At non-vanishing temperatures there are other possible 
quadratically divergent tensor structures. They mainly gi- 
ve different contributions in the infrared which are sub- 
leading as long as the gluon dressing functions are not in- 
frared divergent. The infrared exponents are only weakly 
dependent on the actual choice |22ll5) . Thus we subtract 
the term (f l(if) to remove all spurious divergences. We are 
then left with the intrinsic logarithmic divergences which 
are renormalised in the MOM-scheme as described in ref. 

m 

The ansatze for the vertices © and (|1(J|) are inserted 
into eqs. JHJ and @. After contracting the gluon equation 
with the heat-bath transversal and longitudinal projectors 
© , one arrives at the following three coupled equations for 
the three scalar dressing functions G, Z m and Z$ defined 
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m eqs. 



Q and gj: 



10 



+00 



d z c 



G(q) 



G{k) *' » " J^J (2Ttfk 2 q 2 {k-q) 2 

{A T Z m (k -q) + A L Z {k - q)} , (17) 



1 



Z m (k) 



= Z 3 - - 2 9 



N C Z X T 

N C Z X T 



d 3 q G(q)G(p) 
(2tt) 3 k 2 q 2 p 2 

d 3 q F{q,p,k,Z,G,Z x ) 



R+ 



(2tt) 3 k 2 q 2 p 2 

[M T Z m {q)Z m {p) + M 1 Z Q (q)Z m (p)+ 
M 2 Z (p)Z m (q) + M L Z {q)Z Q {p)} , (18) 

1 7 n 2 N 7T v r d3 i G ^ G ^ 

Z^kj- 23 - 9 NcZlT l^ J ( 27 r)3 t2„2„2 



(2tt) 3 k 2 q 2 p 
d 3 q F{q,p,k,Z,G,Z 1 ) 



P- 



k 2 q 2 p 2 



[N T Z m (q)Z m {p) + N 1 Z (q)Z m {p)- 
N 2 Z (p)Z m (q) + N L Z {q)Z (p)} . 



(19) 



Note that go — w n = lixnT . Herein colour indices have 
been contracted and the tensor S ab has been separated. 
The expressions for the kernel functions At, Al, R, Mt, 
Mi, M 2 M L , P and N T , N x , N 2 , N L can be found in 
appendix A, see eqs. (|32[1 - (I43[1 . A graphical representation 
for these equations is given in fig. ^ 

The MOM-renormalization scheme, applied for van- 
ishing temperature, is used here. This amounts to solv- 
ing subtracted equations, i.e. the respective values at the 
renormalization scale are subtracted on the l.h.s and r.h.s. 
The renormalization conditions are then used to calculate 
the renormalization constants. The introduction of a fi- 
nite temperature does not give rise to any novel diver- 
gences |16j . However, finite, temperature-dependent mod- 
ifications of the renormalization constants occur. Thus we 
impose the same renormalization condition for both the 
transversal gluon propagator and the longitudinal one and 
determine a renormalization constant for each. 

Due to the compactification of space-time, the inte- 
grals become discrete, infinite sums. To solve the equa- 
tions numerically one has to introduce a cut-off A and 
sum only up to this cut-off. The integral equations are 
solved self-consistently. In order to speed up convergence 
a Newton-Raphson iteration is used. As it converges only 
locally we start with results from zero temperature and 
increase the temperature in small steps. The method will 
be explained in more detail elsewhere. 



4 Results 

In order to display the temperature dependence of the 
resulting dressing functions we will show the linear com- 
bination i|12|) of the dressing functions Zq and Z m as well 
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Fig. 2. The ghost dressing function G in the infrared region 
at different temperatures (24 3 momentum grid). 



as their difference 



AZ = Z - Z n 



(20) 



In a first step only the zeroth Matsubara mode is discussed 
here, the other modes will be presented below. The com- 
bination Z in eq. 1)12(1 is the natural one to compare with 
zero temperature results, while AZ measures temperature 
effects on the tensor structure of the gluon propagator. 
The results indicate that G and Z are nearly temperature 
independent, see figs. |21 and |21 The momentum scale is ob- 
tained from the corresponding T=0 calculation |22I ■ Note 
that for the 24 3 grid and the ultraviolet cutoff of 4.7GeV 
corresponds to a spatial volume of approximately 6 fm 3 . 
Unfortunately, it is not possible to extract reliable infrared 
exponents from these torus results. Hence, the changes can 
be either due to changes in the infrared coefficient or in the 
exponent. The ghost dressing would suggest a change in 
the exponent, while the gluon dressing points to a change 
in the coefficient. Nevertheless, the results allow to con- 
clude, that the power laws persist in the infrared regime. 
The changes of the ghost dressing functions with respect 
to temperature are more likely fluctuations due to finite- 
size effects rather than an intrinsic temperature effect. Fig. 
^displays, that AZ is temperature dependent to a rather 
large extent at intermediate momenta. The temperature 
dependences of Z and Z m have opposite signs and AZ 
is of the order of a few percent in the intermediate mo- 
mentum region compared to the sum. Up to about 1 GeV 
the shape of AZ can be explained by just changing the 
infrared coefficient and then connecting continuously to 
the ultraviolet regime, where AZ has to vanish. The data 
presented here are for a cutoff of 4.7 GeV. For other cut- 
offs the results vary slightly, in a range which is expected 
for such low cutoffs. 

For the higher Matsubara modes the dressing functions 
have an interesting property: the combined gluon dressing 
Z and the ghost dressing G seem to depend only on the 
4-momentum (at least approximately), while the higher 
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Fig. 1. Diagrammatic representation of the propagator DSEs in the truncation scheme used here. Wiggly lines denote heat- 
bath transverse gluon propagators, double-dashed lines are heat-bath longitudinal propagators and dashed lines represent ghost 
propagators. Blobs indicate dressed propagators and vertex functions, respectively. 
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Fig. 3. The averaged gluon dressing function Z in the infrared 
region at different temperatures (24 3 momentum grid). 
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Fig. 4. The difference AZ of the gluon dressing functions at 
different temperatures (24 3 momentum grid). 



modes of AZ are not 0(4) invariant, as expected (see 
figs. Ed \ 

Due to the numerical method, we are biased to stay 
in the confining phase. Nevertheless it is surprising that 
we find qualitatively similar solutions also for high tem- 
peratures. These solutions are numerically stable up to 
temperatures allowed by the UV cut-off, 2nT w A. In 
order to draw conclusions on the character of the phase 
transition we need results from other numerical methods 
to support this finding. If so, our results would indicate a 
first order phase transition. Figs. [H] and [5] show that the 
changes of the "wrong-phase" solutions with temperature 
are also smooth and are qualitatively the same as for small 
temperatures. 

The running coupling a(k 2 ), being a renormalization 
group invariant quantity, is of special interest. In Landau- 
gauge at zero temperature, a non-perturbative running 
coupling a(p 2 ) can be defined via the relation: 

a(/)=a(/x 2 )G 2 (p 2 ,/x 2 )Z(p 2 , M 2 ). (21) 

In order to smoothly connect to the results for vanishing 
temperature we take for Z the linear combination l|12|) : 

a(n,plT) :=a(n 2 )G 2 (n,pln 2 ,T)Z(n,pln 2 ,T). (22) 

As expected from the results for the dressing functions, the 
coupling l|22(l shows no temperature dependence beyond 
numerical uncertainties. 

5 CJT Action and Phase Transition 

As we want to investigate the phase transition it is nec- 
essary to compute thermodynamic quantities. In order to 

1 The first Matsubara mode of the ghost dressing (see fig. 01 
deviates slightly from 0(4) invariance. This is due to the torus 
regularization. 
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Fig. 10. Diagrammatic representation of the two-loop contributions to the effective action. 




Fig. 11. ViiPpv, Dg) after inserting the reintegrated DSEs. 




Fig. 5. Different Matsubara modes of the ghost dressing func- 
tion at T = 140 MeV (24 3 momentum grid). 
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Fig. 6. Different Matsubara modes of the averaged gluon 
dressing function ZatT= 140 MeV (24 3 momentum grid) . 



Fig. 7. Different Matsubara modes of the difference AZ of 
the gluon dressing functions at T = 140 MeV (24 3 momentum 
grid). 
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Fig. 9. The averaged gluon dressing function Z at high tem- 
peratures (24 3 momentum grid). 



do so, we need the temperature dependent effective action 
for vanishing fields and full propagators. Since we want to 
calculate the temperature dependence of different quan- 
tities, we need to assure energy conservation. Thus the 
DSEs, which we use to compute the dressing functions, 
have to be the variational equations of the energy func- 
tional, i.e. of the effective action. We take the known 2PI 
effective action by Cornwall, Jackiw and Tomboulis (CJT) 
[H] for fully dressed propagators and bare vertices 2 . First 
we discuss the action for vanishing temperature, given by 
the formula (c.f. ref. 



V(D^,D G ) = V {D^,D G ) + V 2 (D^,D G ) 



(23) 



V (D^,D G ) 



l r 



D, 



d 4 P 

(2^) 4 1 2 



g^\- l -\n(Dl a -\p)D au {p) 
[D%(p)-^D G {j>)-l] 
+ \n{D%{p)^D G {p))} , 



ip)D av (p)- 



(24) 



where we already employed Euclidean space-time conven- 
tions. 

In a truncation scheme with bare vertex functions only 
the diagrams of fig. contribute to the effective action. 
For our truncation to the gluon propagator DSE, the 3- 
gluon vertex (|TU)) is constructed from the dressing func- 
tions of the propagators. However, throughout this inves- 
tigation we neglect variations of the vertex with respect to 



Recently there has been renewed interest in the formalism 
also for higher particle irreducibilities |28|. however, the cor- 
responding treatment is beyond the scope of this paper. Any- 
ways, our truncation scheme only includes 1-loop graphs for 
the propagator DSEs, so the action contains at most two-loop 
diagrams and is thus 2PI. 



the gluon and ghost dressing functions, and thus employ 



(2^ J (2nr ipMp - q) 
k»D G (q)D G (p) -j^J T^fl/ tJ^I { 
r ° pa (P> -P - 1i q) D af}{q)D pa {p + q) 
r° al/ (-q,p + q,-p)D flu (p)} . (25) 

The variation of the action l)23|) with respect to the propa- 
gators reproduces the Dyson-Schwinger equations for the 
gluon and ghost propagators © and (jSJ), however, with 
bare 3-gluon vertex. Multiplying the gluon DSE with the 
full gluon propagator and the ghost DSE with the full 
ghost propagator and performing the trace in colour and 
Lorentz indices and the integral in momentum space, we 
recover the two-loop diagrams of the CJT action. This 
results in: 



V 2 (D^,D G ) 



d' P r 1 r 

(2^) 4 \ 6 



f, a {p)D au {p)- 

M + \ [D%{p)- l D G {p) 1] } , (26) 

see fig. ^2 for a graphical representation. Inserting (|26[l 
and the explicit expression for the propagators in depen- 
dence on the dressing functions in l|23(l one obtains: 



V(Z,G) = N C J ^2vr) 4 { 



{Z 3 Z(p 2 ) 
2 



-In (Z 3 Z(p 2 )) 



- (z 3 G(p 2 ) In (z 3 G(p 2 ))}. (27) 

For solving the corresponding DSEs numerically, we in- 
troduce a momentum cut-off A. Thus the renormalization 
constants Z 3 for the gluon wave function and Z 3 for the 
ghost propagator not only depend on the renormalization 
scale fj, but also on the cut-off 3 . Defining x :— p 2 one has: 



V(Z,G) 



N c 
(4tt) 2 
3 

2 
2 



dxx 



{(Z 3 (»,A)Z(x,ii) 



1) 



3 



]n(Z 3 (n,A)Z(x,n)) 
Z 3 ( f i,A)G(x, t x)-l 



hi 



(z 3 0M)G(z, M )) }• 



(28) 



This quantity diverges like yl 4 , if the renormalization con- 
stants are only known to a finite accuracy. Furthermore, 
the three-gluon vertex is dressed by construction, which 
is not accounted for by the above action. In fact the fol- 
lowing fit 



fix) 



ax 



(29) 



3 The propagators are cut-off independent within numerical 
accuracy as soon as A > 15 GeV. 
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to the cut-off dependence of the action yields: 

a = -3.2 x 1(T 3 

c = (9.8 ± 1.3) x l(T 10 GeV 4 



(30) 



The result should be interpreted as vanishing of the effec- 
tive action as the first term is of purely numerical origin 
and the second term is much too small to correspond to a 
scale of the system, « 5.6 MeV. 

At finite temperatures, the integral measure in eq. H27J1 
has to be replaced by that for the space-time torus, see eq. 
(T7J|. Furthermore, the finite temperature gluon propagator 
with its additional tensor structure 10 has to be inserted. 
Finally we get for the temperature dependent 1-particle 
irreducible CJT- Action: 



V(T) = N C T 



+ 00 

£ 

n= — oo 



d 3 q 

(27T) 



3 { 



Ki |p|) - 1) - In (Z 3m Z m (to n , \p\)) + 
(Z 30 Z K, |p|) - 1) - i In (Z 3Q Z (uj n , |p|)) - 

[z 5 G{u n , |p|) - l) + In (z a G(u n , |p|)) } . (31) 

This effective action is temperature independent within 
numerical uncertainties up to T = 250 MeV. For higher 
temperatures large variations are seen. Whether this is 
connected to the phase transition remains to be seen. 

6 Conclusions 

In this paper we have extended a truncation scheme for 
the Dyson-Schwinger equations of Yang-Mills theories in 
Landau-gauge to non-vanishing temperatures using the 
imaginary time formalism for quantum field theories in 
thermal equilibrium. We have obtained numerical solu- 
tions for the gluon and ghost propagators. Furthermore 
we have derived an expression for the CJT action of the 
interacting Yang-Mills theory, depending on dressed prop- 
agators. 

The results can be summarised as follows: Tempera- 
ture dependences of the ghost and gluon propagators are 
rather weak and the power-law behaviour in the infrared, 
already observed at vanishing temperature, persists. The 
corresponding CJT action vanishes (within numerical un- 
certainties) and does not possess any significant temper- 
ature dependence up to T — 250 MeV. Thus gluons and 
ghost neither build up any pressure nor give raise to finite 
entropy of the system. An improved numerical method, es- 
pecially one suitable for infinite volume, is, of course, desir- 
able to study thermodynamic quantities below the phase 
transition. Corresponding investigations are currently per- 
formed. 
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Appendix 

A Integral kernels for the DSEs 

The complete algebraic expressions for the Lorentz struc- 
tures that occur in the DSEs for the different products of 
dressing functions, are given here. First, there are the two 
kernels of the ghost equation l|17fl : 



A T (k,q) = - 



k 2 q 2 -{k-q) 
(k-q) 2 



(32) 



A,( fc , g )=- fc2 ?:- ( y + fc x- (fc ; g)2 - (33) 



(k - qf 



(k-q) 



Secondly, the kernel of the ghost loop and the four kernels 
of the gluon loop for the heat-bath transversal part of the 
gluon equation {TSJ: 



( g 2 fc 2 - (fc ■ qf) 
R ( k >l) = p 

M T (k, q) = ^l;^^ 2 ((fc • qf + k 2 q 2 - 



(34) 





2p 2 (fc 2 + q 2 )) 




(35) 


M 1 (k,q) 


(<?ofc ■ q - 


- k q 2 ) 2 ((fc • p) 2 + fc 2 p 2 ) 


(36) 




k 2 q 2 p 2 q 2 


M 2 (k,q) 


(pofc • P 


-k p 2 ) 2 ((k-q) 2 + k 2 q 2 ) 


(37) 




k 2 q 2 p 2 p 2 


M L {k,q) 


{q 2 k 2 - 
- 2 vy 


{k-q) 2 ) (q 2 p 2 -qopmpf 


(38) 




k 2 q 2 p 2 k 2 p 2 



And third the corresponding kernels for the heat-bath lon- 
gitudinal part (|P31) : 



P(k,q) 

N T (k,q) 

N 1 (k,q) = -2 

N 2 (k,q) 
N L (k,q) 



(q k 2 - fc fc • q) 



k 2 k 2 

(k k ■ q - q k 2 ) 2 ((q ■ p) 2 + q 2 p 2 



k 2 q 2 p 2 k 2 




o (q 2 k 2 -(k-q) 2 ) (k 2 q 2 


- k q kq) 


k 2 q 2 p 2 k 2 p 2 




Jq 2 k 2 -(k-q) 2 ) (k 2 p 2 


- k Q p Q kp) 2 



2k 2 q 2 p 2 k 2 q 2 p 2 
(p • q)k 2 (p q 2 + qop 2 ) 
(fc • q)p 2 {k Q q 2 - q k 2 )] 



k 2 q 2 p 2 k 2 q 2 

[(p-fc)q 2 (k p 2 +p Q k 2 ) 



(39) 
(40) 
(41) 
(42) 



(43) 
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Fig. 12. The ghost dressing function G at T 
different momentum grids. 



f 40 MeV on Fig. 14. The difference AZ of the gluon dressing functions at 
T = 140 MeV on different momentum grids. 




k 3 [GeV] 



Fig. 13. The averaged gluon dressing function Z at T = 140 
MeV on different momentum grids. 



B Finite size effects 

As expected, finite-size effects cannot be neglected com- 
pletely but are rather small and change the results only 
in a quantitative way, as can be seen from the figures l"H?l 
- ^] The deviations tend to increase for higher tempera- 
tures, since the ultraviolet cutoff is rather low (4.7 GeV). 
Finite size effects do not alter our conclusions. Especially, 
the temperature estimate of eq. iJSJ is fairly independent 
of the finite-momentum lattice size. 
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